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Abstract 

The motion of a rigid body in a quadratic potential is an important example of an integrable 
Hamiltonian system on a dual to a semidirect product Lie algebra so(n) x Symm(n). We give a 
Lagrangian derivation of the corresponding equations of motion, and introduce a discrete time 
analog of this system. The construction is based on the discrete time Lagrangian mechanics 
on Lie groups, accompanied with the discrete time Lagrangian reduction. The resulting multi- 
valued map (correspondence) on the dual to so(n) k Symm(n) is Poisson with respect to the 
Lie-Poisson bracket, and is also completely integrable. We find a Lax representation based on 
matrix factorisations, in the spirit of Veselov-Moser. 
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1 Introduction 



The rigid body dynamics are rich with problems interesting from the mathematical point of view, 
in particular, with integrable problems. Certainly, the most famous ones are the three integrable 
cases, named after Euler, Lagrange, and Kovalevskaya, of the rotation of a heavy rigid body around 
a fixed point in a homogeneous gravity field. They were discovered in the 18th and the 19th century, 
and can be called classical. However, the list of intergrable problems of the rigid body dynamics is 
by far not exhausted by these ones. 

In the present paper we turn our attention to the rotation of a rigid body around its fixed 
center of mass in an arbitrary quadratic potential. The integrability of this problem is a much 
more recent observation due to Rcyman [R] and Bogoyavlensky [B]. (However, some particular 
case of this result was given already by Brun [Br] ; the equations of motion in this case are identical 
with those describing the integrable case of the motion of a rigid body in an ideal fluid, due to 
Clebsch [C]). 

A problem we solve in the present paper, is a construction of an integrable discretizatAon of 
the above mentioned integrable mechanical system. Our construction is close in spirit to the work 
by Moser and Veselov [MV], where they used an approach based on the discrete time Lagrangian 
mechanics. Several important integrable mechanical problems were discretized in [MV], including 
the Euler top and the Neumann system. The general theory of discrete time Lagrangian mechanics 
on Lie groups was developed further in [BSl], where the list of Moser and Veselov was extended 
by a discrete time Lagrange top. A further development of the general theory was undertaken in 
[BS2], where general discrete time Euler-Poisson equations on semidirect product Lie algebras were 
obtained as a result of a redution procedure, applied to discrete time Lagrangian systems on Lie 
groups. A discretization of a top in a quadratic potential, achieved in the present paper, serves as 
a spectacular illustration to the abstract constructions in [BS2]. 

We recall the general theory of the Lagrangian reduction in the continuous time and in the 
discrete time contexts, respectively, in Sect. 2 and 3. Further, we give in Sect. 4 a Lagrangian 
derivation of the equations of motion of a rigid body in a quadratic potential. In Sect. 5 we 
introduce the discrete time analogs of these constructions. Finally, in Sect. 6 the specialization of 
these results for the Clebsch case is given. Conclusions are contained in Sect. 7. 

2 Lagrangian mechanics and Lagrangian reduction on TG 

Recall that a continuos time Lagrangian system is defined by a smooth function J-i{g, g) : TG i— > M 
on the tangent bundle of a smooth manifold G. The function L is called the Lagrange function. 
We will be dealing here only with the case when G carries an additional structure of a Lie grioup. 
For an arbitrary function g(t) : [ioj^i] '-^ G one can consider the action functional 



J to 

A standard argument shows that the functions g{t) yielding extrema of this functional (in the class 
of variations preserving g{tQ) and g{ti)), satisfy with necessity the Euler-Lagrange equations. In 




(2.1) 
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local coordinates {g^} on G they read: 

d fdL\ dL 



dt ydg' ) dg' ' ^^'^^ 

The action functional S is independent of the choice of local coordinates, and thus the Euler- 
Lagrange equations are actually coordinate independent as well. For a coordinate-free description 
in the language of differential geometry, see 0], |MR]. 
Introducing the quantities 

n = VgL € T*G , (2.3) 

one defines the Legendre transformation: 

{g,g)eTG ^ {g,U)GT*G. (2.4) 

If it is invertible, i.e. if g can be expressed through {g, 11), then the Legendre transformation of the 
Euler-Lagrange equations ( |2.2| ) yield a Hamiltonian system on T*G with respect to the standard 
symplectic structure on T*G and with the Hamilton function 

H{g,n) = {U,g)-L{g,g) , (2.5) 

(where, of course, g has to be expressed through {g,lV)). 

When working with the tangent bundle of a Lie group, it is convenient to trivialize it, translating 
all vectors to the group unit by left or right multiplication. We consider only the left trivialization 
here: 

{g,g)eTG ^ {g,n)eGxg, (2.6) 

where 

n = Lg-iJ g = Lg^n . (2.7) 

Denote the Lagrange function pushed through the above map, by L*^')((7, 17) : G x g i— > M, so that 

L^'Hg,n)=L{g,g) , (2.8) 



The trivialization (|2.6| ) of the tangent bundle TG induces the following trivialization of the 
cotangent bundle T*G: 

{g,U)eT*G ^ {g,M)eGxQ* , (2.9) 

where 

M = L*gU ^ U = L*-iM. (2.10) 
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We now consider the Lagrangian reduction procedure, in the case when the Lagrange function is 
symmetric with respect to the left action of a certain subgroup of G, namely an isotropy subgroup 
of some element in the representation space of G. So, the next ingredient of our construction is a 
representation ^ : G xV i-^V oi a Lie group G in a linear space V; we denote it by 

$(5) -v for g e G , V eV . 

We denote also by the corresponding representation of the Lie algebra gmV: 



for ^Gfl, veV. (2.11) 

e=0 



The map (f>* : q x V* >->■ V* defined by 

{<t>*{0-y,v) = {y,m-v) yveV , yeV* , ^€Q, (2.12) 

is an anti-representation of the Lie algebra g in V* . We shall use also the bilinear operation 
o : V* X V i-^ g* introduced in [HMR, CHMR] and defined as follows: let v G V , y G V* , then 

{yov,0 = -{y,m-v) v^Gfl. (2.13) 

(Notice that the pairings on the left-hand side and on the right-hand side of the latter equation 
are defined on different spaces). 

Fix an element a eV, and consider the isotropy subgroup G^"! of a, i.e. 

= {h : $(/i) ■a = a}cG . (2.14) 

Suppose that the Lagrange function Ij(g,g) is invariant under the action of G^"] on TG induced by 
left translations on G: 

L{hg,Lh,g) = Ug,g) , ^ G G^ . (2.15) 

The corresponding invariance property of 'L^^\g,Q) is expressed as: 

L^^\hg,n) ='L^^\g,n) , /i G G^ . (2.16) 

We want to reduce the Euler-Lagrange equations with respect to this left action. As a section 
(G X fl)/Gl"] we choose the set g x Oa, where Oa is the orbit of a under the action $: 

Oa = mg)-a, geG}cV . (2.17) 

The reduction map is 

{g,n) eGxg ^ (^^, P) G X 0„ , where P = ^ig~^) ■ a , (2.18) 
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so that the reduced Lagrange function C^^'^ : g x Oa ^ M is defined as 

£(')(17,P) = L(')(5,0) , where P = ^[g-^)-a. (2.19) 
The reduced Lagrangian C^^\n,P) is weh defined, because from 

P = ^g^') ■ a = <^{g^^) ■ a 
there fohows <I>(<^25]" ) • a = a, so that 525r^ e G''"^' and L«(5i,17) = L«(52,5^). 

Theorem 1 [HMR, CHMR] 

a) Under the left trivialization {g,g) i— > (g,^) and the subsequent reduction (g,^) the 
Euler-Lagrange equations (2.2) become the following Euler— Poincare equations; 

M = ad*n-M + VpC^^^ o P , 

(2.20) 

p = -<Pin)-p, 

where 

M = L*n = Vn/:^'^ . (2.21) 

b) If the "Legendre transformation" 

(n, P) G g X Oa ^ (M, P) G g* X , (2.22) 
is invertible, then ( 2.20 ) is a Hamiltonian system of q* x Oa with the Hamilton function 

H{M, P) = (M, Q) - {n,P) , 
with respect to the Poisson bracket given by 

{Pi, P2} = (M, [VmFi, VmF2] ) + (VpPi, (^(VAfP2) • P) - (VpP2, <^(VMi^i) -P) 

(2.23) 

for two arbitrary functions Pi 2(M, P) : g* x Oq ^ M. 

Remark 1. The formula ( |2.23| ) defines a Poisson bracket not only on g* x Oq, but on all of 
g* X Rewriting this formula as 

{Pi, P2} = (M, [VmPi, VmF2] ) + {P, (I^*{^mF2) ■ VpPi - r (Va/Pi) • VpP2 ) 

(2.24) 

one immediately identifies this bracket with the Lie-Poisson bracket of the semiproduct Lie algebra 
Q t< V* corresponding to the representation —(p* of g in V*. 

Remark 2. In an important particular case of constructions of this section, the vector space 
is chosen as the Lie algebra of our basic Lie group: V = q, the group representation is the adjoint 
one: $(5) • v = Ad g ■ v, so that (/>(^) ■ v = ad ■ v = [£,,v], and the bilinear operation o is nothing 
but the coadjoint action of g on g*: yov = ad* v ■ y. This is the framework, e.g., for the heavy top 
mechanics. 
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3 Lagrangian mechanics and Lagrangian reduction on G x G 



We now turn to the discrete time analog of these constructions, Our presentation of the general 
discrete time Lagrangian mechanics is an adaptation of the Moser-Veselov construction [V] , [MV] 
for the case when the basic manifold is a Lie group. The presentation of the discrete time La- 
grangian reduction follows [BSl], [BS2]. Almost all constructions and results of the continuous 
time Lagrangian mechanics have their discrete time analogs. The only exception is the existence 



of the "energy" integral (|2.5D . 

Let Li{g,g) : G x G i— > M be a smooth function, called the (discrete time) Lagrange function. 
For an arbitrary sequence {gk € G, k = ko, ko + 1, . . . ,ki} one can consider the action functional 

§=Y,U9k,9k+i)- (3.1) 

k=ko 

Obviously, the sequences {gk} delivering extrema of this functional (in the class of variations 
preserving g^^ and gki), satisfy with necessity the discrete Euler-Lagrange equations: |^ 

ViL(5-fc, fiffc+i) + V2L{gk-u9k) = . (3.2) 

Here ViL((7, g) (V2L((7, g)) denotes the gradient of L((7, g) with respect to the first argument g (resp. 
the second argument g). So, in our case, when G is a Lie group and not just a general smooth 



manifold, the equation (3.2) is written in a coordinate free form, using the intrinsic notions of the 
Lie theory. As pointed out above, an invariant formulation of the Euler-Lagrange equations in the 
continuous time case is more sophisticated. This seems to underline the fundamental character of 
discrete Euler-Lagrange equations. 

The equation (|3.2| ) is an implicit equation for gk+i- In general, it has more than one solution, 
and therefore defines a correspondence (multi-valued map) {gk-i,9k) ^ i9k,9k+i)- To discuss 
symplectic properties of this correspondence, one defines: 

Uk = S72Ugk^u9k) eT*^G . (3.3) 

Then may be rewritten as the following system: 

Ilfc = -Vih{gk,gk+i) 

(3.4) 

rife+i = V2L(gfc,5fc+i) 



^For the notations from the Lie groups theory used in this and subsequent sections see, e.g., [BSl]. In particular, 
for an arbitrary smooth function / : G i-^ R its right Lie derivative d'f and left Lie derivative df are functions from 
G into 0* defined via the formulas 



V77G. 



and the gradient Vf{g) G TgG is defined as 

^fig) = R;-^df{g)^L;-^d'f{g). 
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This system defines a (multivalued) map (^fc, Ilfc) i— > (^fc+i, Ilfc+i) of T*G into itself. More precisely, 
the first equation in ( |3.4| ) is an implicit equation for gu+i, while the second one allows for the explicit 
and unique calculation of Ilfc+i, knowing and gu+i- As demonstrated in 0, | MV ], this map 
T*G ^ T*G is symplectic with respect to the standard symplectic structure on T*G. 

The tangent bundle TG does not appear in the discrete time context at all. On the contrary, 
the cotangent bundle T*G still plays an important role in the discrete time theory, as the phase 
space with the canonical invariant symplectic structure. The left trivialization of T*G is same as 
in the continuous time case: 



where 



Consider also the map 



where 



{gk,Iik) G T*G ^ igk,Mk) G G x g* , (3.5) 



Mk = L;Uk ^ Uk = L*-,Mk. (3.6) 



{gk,9k+i) eGxG ^ {gk, Wk)£GxG, (3.7) 



Wk = ^gk+i ^ gk+i = gkWk . (3.8) 

Denote the Lagrange function pushed through ( p. 7] ) by 

h^'\gk,Wk)=Ugk,gk+i) ■ (3.9) 

Suppose that the Lagrange function I^{g,g) is invariant under the action of Gt"' on G x G 
induced by left translations on G: 

Uhg,hg) = Ug,g) , /i G G^ . (3.10) 

The corresponding invariance property of 'L^^^g, W) is expressed as: 

lW(/i5,VF) =lW(5,VF) , /igGI"]. (3.11) 

We want to reduce the Euler-Lagrange equations with respect to this left action. As a section 
(G X G)/GM we choose the set G x O^. The reduction map is 

{g,W)eGxG ^ {W,P)eGxOa, where P = ^{g-^) ■ a , (3.12) 

so that the reduced Lagrange function A^'^ : G x 1^ is defined as 

AW(H^,P) =l(')(5,VF) , where P = ^{g~^) ■ a . (3.13) 
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Theorem 2 [BSl], [BS2] 

a) Under the left trivialization {g,g) i— > {g,W) and the subsequent reduction {g,W) (W,P), 
the Euler-Lagrange equations ( |3.2D become the following discrete Euler— Poincare equations; 

Ad* W^:^ ■ Mk+i = Mfc + VpA«(H^fc, Pfc) o Pk , 

(3.14) 

Pk+i = '^{W^^) ■ Pk , 
where 

Mk = d{vk^'\Wk-i.Pk-i) G 5* . (3.15) 

b) If the "Legendre transformation" 

{Wk-i,Pk-i) £GxOa^ {Mk,Pk) € B* X 0„ , (3.16) 

where Pk = ^{Wjr\) ■ Pk-i, is invertible, then ( |3.14| ) define a map {Mk,Pk) i— > (^fc+i, -Pfc+i) of 
0* X Oa which is Poisson with respect to the Poisson bracket ( p.23| ). 

The relation between the continuous time and the discrete time equations is established, if we 

set 

ak = g, gk+i = g + eg + 0{e^) , l.{gk,gk+i) = eM9,g) + 0{e^); 
Pk = P, Wk = l + en + 0{£^) , A«(W^fe,Pfc) =e£«(l^,P) + 0(e2) . 

4 A rigid body in a quadratic potential 

The basic Lie group relevant for our main example is 

G = SO(n) , so that q = so(n) 
(the "physical" rigid body corresponds to n = 3). The scalar product on g is defined as 

(C,??) = -^tr(^r?) , C,?7G0. 

This scalar product is used also to identify q* with g, so that the previous formula can be considered 
also a pairing betwenn the elements ^ G and r/ G g*. 

The group G is a natural configuration space for problems related to the rotation of a rigid 
body. Indeed, if E{t) = ^ei(t), . . . , e„(t)^ stands for the time evolution of a certain orthonormal 
frame firmly attached to the rigid body (so that all Ck G M"), then 

E{t) = g~\t)E{0) ^ ek{t) = g'\t)ek{0) {I < k < n) , 
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with some g{t) G G. The Lagrange function of an arbitrary rigid body rotating abot a fixed point 
G M" in a field with a quadratic potential 



1 " 

fix) = 9 XI '^ij^i^j ' (4-1) 

is equal to (cf. [B], [RSTS]): 



2 

2j = l 



Ug, g) = -ltv{njn) + ^ iv{gJg^A) , (4.2) 



where 



= g g = g g is the angular velocity of the rigid body in the body frame E{t); 



• J is a symmetric matrix (tensor of inertia of the rigid body); choosing the frame E{0) properly, 
we can assure this matrix to be diagonal, J = diag(Ji, . . . , J„), which will be supposed from 
now on; 

• A = (a.jj)'"^^^ is a symmetric matrix of coefficients of the quadratic form ip{x). 

To include this Lagrange function in the framework of Sect. 2, we make the following identifications: 

• V = Symm(n), the linear space of all n x n symmetric matrices; we identify V* with V via 
the following scalar product on V: 

{VI,V2) = ^tT{viV2) , Vl,V2eV. 

• The representation $ of G in F is defined as 

^(g) ■ ^ = g'vg~^ = g^g^ fo^ g g G, v g v . 

• Therefore the representation of g in F is given by 

^iO ■v = [C,v] for ^eg, veV , 
while the anti-representation of g in V* is given by 

<l>*{0-y = -[ty] for e eg, ye T^*, 

• Finally, the bilinear operation o : F* x g* is given by 

y ov = —[y,v] ioi y E V* , V E V . 
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Denoting now P = Ag = ^{g ^) ■ A, we represent ( [4 .21) in the form 

L(ff, g) = £(')(0, P) = tv{njn) + ^ tr( JP) = ti{njn) + (J, P) , (4.3) 

which is manifestly invariant under the left action of the isotropy subgroup G^-^^ . Now Theorem 1 
is applicable, which delivers the following equations of motion: 

M = [M, n] + [p, J] , 

(4.4) 

p = [p,n] , 

where 

M = VnC^^^ = jn + nj ^ Mjk = {Jj + Jk)^jk ■ (4.5) 



According to the general theory, the system ( |4.4| ) is Hamiltonian on the dual of the semidirect 
product Lie algebra q t< V*, with the Hamilton function 

1 1 1 " 

HiM, p) = -{M, n) -{J,P) = ^Y1 -rih- - 9 E 'J^P^'^ ■ (4-6) 
^ ^ j<k -^^ + ^ k=l 

Generic orbits in this Poisson phase space have dimension n'^ — n (the dimension of g x y* is equal 
to n(n — l)/2 + n(n + l)/2 = n^; the n spectral invariants of P E y are Casimir functions of the 
Poisson bracket). In particular, for the "physical" case n = 3 the dimension of the generic orbit 
is equal to 6. Therefore the number of independent involutive integrals of motion necessary for 
complete integrability is equal to n(n — l)/2 (equal to 3 for n = 3). 

A key observation of [R], [B] consists in the following Lax representation of the above system: 

L{X) = [L{X),B{X)], 

where 

L(A) = P + AM + B{X) = n + XJ . 

(The key point in the straightforward verification of this statement is the identity 

[M, J] + [j2,f)] = 0, 

which follows directly from (4.5)). The spectral invariants of the matrix L{X) provide us with the 
necessary number of independent integrals of motion [B], and their involutivity follows from the 
general r-matrix theory (cf. [R], [RSTS]). 
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5 A discrete time analog of a top in a quadratic potential 

To find a discrete analog of the Lagrange function ( |4.2D , we rewrite the latter once more as 

L{g, g) = ^ tT{gjf ) + \ irigJg^A) . (5.1) 

Let us introduce the following discrete analog: 

M9k,9k+i) = ^tr(^(c/fe+i -gk)J{gk+i - 5fc)'^) + | ^(5-^+1 J^JA) . (5.2) 

The powers of e are introduced in a way assuring the correct asymptotics of the discrete time 
Lagrange function, namely ^{gk, Ok+i) ~ ^^{g,g), Qk = 9 and gk+i ~ 9 + ^9 (see the end of Sect. 
3). Up to a constant, the function ( ^.2| ) may be rewritten as 

M9k,9k+i) = -^i^igk+iJgl) + ^^^{9k+iJ9k A) ■ (5.3) 

This is representable also in terms of Wk = gjgk+i € G and Pk = g'k^gk £ Oa C V: 

H9k,9k+i) = A^'HWk,Pk) = --^tiiWkJ) + ^tT{WkJPk) . (5.4) 

(Recall that in the continuous limit one has Wk ~ / + eft, so that ~ / — eQ). 

Theorem 3 The discrete time Euler-Lagrange equations for the Lagrange function ( ^.4[ ) are equiv- 
alent to the following system: 

Mk = \{WkJ - JW^) - ^{PkWkJ - JW^Pk) 

' Mk+i = -{JWk-W^J)-UjPkWk-W^PkJ) , ^^-^^ 
^ Pk+i = W^PkWk . 

The multi-valued map (correspondence) {Mk, Pk) ^ {Mk+i, Pk+i) described by ( |5.5| ) is Poisson with 
respect to the Lie-Poisson bracket of the semidirect product Lie algebra q t< V* , where q = so{n), 
V* = Symm(n), and the representation —(p* of g in V* is defined as —(p*i(,) ■ y = [C^y] for C ^ 9> 
yeV*. 

Proof. We are in a position to apply Theorem 2. To this end we first calculate 

VpA^'\Wk,Pk) = ^{WkJ + JWl), 
so that the equations of motion read: 

' WkMk+iW^ = Mk + \ [Pk, WkJ + JW;^] 

(5.6) 

. Pk+i = W^PkWk . 
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These equations are, obviously, a discrete time approximation of ([4.4|). They have to be supple- 
mented by a discrete version of (|4.5|), which reads: 



Mfc+i = dl^A^'HWk, Pk) = -{JWk - W'^J) - ^{JPkWk - W;^PkJ) . (5.7) 



Plugging this into the first equation of the system (5.6), we put this system into the form (|57 

The definition of the above correspondence ( ^.5| ) crucially depends on the solvability of the first 
equation in ( |5.5D for € G. The best approach to this problem, as well as to the integrability of 
the correspondence, is through the matrix factorizations. The following argument is a generalization 
of the Moser-Veselov approach [MV] to the discrete time motion of the free n-dimensional rigid 
body. 

Let us introduce the following matrices, depending on the spectral parameter A: 

Uk{X) = [l - ^ Pk)Wk - eXJ , (5.8) 



2 



SO that 

Uj{-X) = W^[l-'-Pk)+ eXJ . (5.9) 
Denote also 

Lfc(A) = (/ - ^ Pfc) ' - e^XMk - e^AV^ = I - (p,, + XMk + AV^) + ^ P| . 

^ ^ ^ ^ ^ ^ (5.10) 

Then a direct calculation allows us to verify the following statement. 

Theorem 4 The equations of motion ( |5.5| ) are equivalent to the following matrix factorizations: 

Lk{X)=UkiX)Uji-X) , 

(5.11) 

Lk+i{X)=U^{-X)Uk{X) . 
In particular, the matrix Lk{X) remains isospectral in the discrete time evolution described by the 



equations (5.5). 



More precisely, the first equation in (|5.5|) is equivalent to the first equation in ( ^.11 ) (under the 



assumption that Wk € G = SO(n)), while the second and the third equations in (|5.5D are equivalent 



to the second equation in ( 5.11 ). So, the problem of solvability of the first equation in (|5.5| ) for Wk 
is equivalent to the matrix factorization problem expressed by the first equation in ( 5.11| ). For a 
treatment of a closely related factorization problem, we refer the reader to [MV]. 

From Theorem 4 there follows also the complete integrability of our discrete time Lagrangian 
map. Notice that its integrals of motion do not coincide with the integrals of motion of the 
continuous-time problem (with the only exception of the free rigid body motion considered in 
[MV]). To be more concrete, the integrals of motion of our map are obtained from the integrals of 
the continuous time problem by replacing P through P — je^P^. As for the actual integration of 
our map in terms of theta-functions, we leave it as an open problem for want of a better occasion 
(cf. [MV], [B], [RSTS]). 
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6 A particular case: the Clebsch problem 

An important particular case of the rigid body in a quadratic potential appears when A = aa^ 
with some a G M". For example, this is the case when the quadratic potential (4.1) represents 
the quadratic terms in an expansion of a potential of a single point mass; it is supposed that the 
distance from the rigid body to this point mass is much larger than the size of the body itself, and 
the ratio of this two length scales is the small parameter of the above mentioned expansion. In this 
case the vector a points from the point mass to the fixed center of mass of the rigid body. Then 

p = g^a (6.1) 

represents the same vector in the frame firmly attached to the rigid body. We have: 

P = pp^ , (6.2) 

i.e. the orbit Oa in V consists of rank 1 matrices. This case could be considered independently, 
along the lines of this paper. The relevant representation space V would be then M", and the 
representation $ of G in F would be defined as ^{g) ■ v = gv. However, we prefere to simply use 
the results already obtained, replacing A with aa^, and P with pp^. So, the Lagrange function 
([4.3| ) takes the form 

L{g,g) = C'^'\n,p) = -^tv{njn) + ^{p,Jp) , (6.3) 

where now (•, •) denote the standard scalar product in M". The equations of motion (4.4) become 

M = [M,n] +pA{Jp) , 

(6.4) 



p 



-VLp 



where the notation pAq = pq^ — qp"^ for p,q G M" is used. This system is Hamiltonian on the dual 
of the semidirect product Lie algebra e(n) = so(n) k M", with the Hamilton function 



1 



1 



H{M,p) = {M,n)--{p,Jp) 



1 



2 Jj + Jk 



1 



(6.5) 



j<k •' k=l 

The integrable discretization of this system is given by the discrete time Lagrange function 



h{gk,gk+i) = A^^\Wk,Pk) 



-- tT{WkJ) + I {pk+i, Jpk) 



--iT{WkJ) + I {pk,WkJpk) , 



(6.6) 



where, as usual, 
tion read: 



Qkdk+i S G and pk = gja G Oa- The equations of motion of this discretiza- 
1 



Mk 



Mk 



-{WkJ 



■Pk A {Jpk 



+1) 



+1 



-{JWk 



W^J) + ^Pk+iAiJpk) 



(6.7) 



Pk+i = W^pk 
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The Lax matrix and the Lax representation of this map are obtained from ( |5.1[l| ) and Theorem 5 
by replacing Pk through PkPk- 

It remains to be noticed that originally the integrable Hamiltonian (|6.5D was found (in the 
"physical" case n = 3) by Clebsch [C] in another setting, namely in the problem of a motion of 
a rigid body in an ideal fluid. (An n-dimensional generalization is due to Perelomov [P]). In this 
setting the system appears as a result of reduction of a Lagrangian on the group E(n), left-invariant 
under the action of a whole group, rather than a Lagrangian on SO(n), left -invariant under the 
action of an isotropy subgroup of a. Nevertheless, these two different settings lead to formally 
identical results. 



7 Conclusion 

The model introduced in the present paper serves as a further important example of the completely 
integrable Lagrangian systems with a discrete time a la Moser-Veselov. Actually, this is the second 
example (after [BSl]) where the version of the discrete Lagrangian reduction is essential, leading 
to systems on duals to semidirect product Lie algebras. This new application was made possible 
due to the theoretical development in [BS2]. Probably, the type of Lagrangians introduced in these 
papers is able to produce further interesting examples, important for applications. At this point, 
I would like to express my gratitude to A.Bobenko, the collaboration with whom in [BSl], [BS2] 
was crucial also for this work. 

Generally, we consider the discrete time Lagrangian mechanics as an important source of sym- 
plectic and, more general, Poisson maps. From some points of view the variational (Lagrangian) 
structure is even more fundamental and important than the Poisson (Hamiltonian) one (cf. pMIl| , 



I MPS I , where a similar viewpoint is represented) . A very intriguing and still not completely under- 
stood point is a capability of the discrete Lagrangian approach to produce completely integrable 
systems. It would be highly desirable to continue the search for integrable Lagrangian discretiza- 
tions of the known integrable systems. Hopefully, the list of such discretizations established in 
[V],[MV], [BSl], and the present paper, will be further extended. Also the generalizations to the 
infinite dimensional case, e.g. to discretization of ideal compressible fluids motion (see [HMR]), are 
highly desirable. 
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